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In  recent  years  there  has  been  considerable  experimental  and 
some  theoretical^’^>interest  in  the  spatial-temporal  growth  of  instabilities 
arising  when  a  beam  is  injected  into  a  plasma.  There  does  not  appear,  how¬ 
ever,  to  be  any  theoretical  treatment  of  the  modification  of  the  growth  of 
an  initial  perturbation  due  to  the  presence  of  a  spatial  boundary.  This 
modification  will  be  illustrated  for  the  case  of  a  monoenergetic  electron 
beam  injected  into  a  cold  plasma  with  an  infinitely  massive  ion  background. 

The  treatment  is  based  on  the  linearized  hydrodynamic  equations 
and  Maxwell's  equations  for  the  normalized  densities,  the  drift  velocities. 
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Eq.  (2)  contains  a  phenomenological  collision  term,  representing  collisions 
between  the  plasma  electrons  and  the  ion  background.  The  remaining  Maxwell's 
equation: 

e,x  +  «P2|*P  +  “U^B  *  °  (6) 

is  treated  as  a  condition  to  be  satisfied  at  t  *  0,  whence  it  is  satisfied  for 
all  later  times. 

Eqs.  (1-5),  which  are  of  the  following  form  (repeated  indices  are 


summed): 
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are  to  be  solved  for  x  >  0  and  t  >  0  subject  to  conditions:  f^(0,t)  *  g^(t) 
and  f^(x,0)  =  h^x).  These  conditions  are  not  ai*  independent  since  Eqs.  (2) 
and  (5)  must  be  satisfied  at  x  =  0,  and  equation  (6)  must  be  satisfied  at 
t  *  0.  The  first  two  of  these  restrictions  are  due  to  the  fact  that  char¬ 
acteristics  of  Eqs.  (1-5)  lie  on  the  boundary  x  *  0. 

The  double  Laplace  transformation  solution  to  Eqs.  (1-5)  is: 

fi(x,t)  *  (2rri)  2J*dpdq  ept  +  qx  Fi(p,q)  (8) 

^(p.q)  “  CajkHk<q>  +  bjkGk<p>3  dji(p.q)/D(?»q) 

and  are  the  one  dimensional  Laplace  transforms  of  h^(x)  and  g^(t). 

D  is  the  determinant  of  the  matrix  a..p  +  b..q  +  c..,  and  d. .  is  the  cofactor. 

ij  iJ  ij 

(D  *  0  is  the  dispersion  relation,  apart  from  a  multiplicative  factor  p,  for 
the  infinite  case.)  The  integration  paths  i  i  (8)  are  to  the  right  of  all 
singularities  of  the  integrand. 

If  the  transforms  contain  only  simple  poles  at  q  *  Of.,  i  *  1..., 
the  q  integration  in  Eq.  (8)  may  be  performed  to  gives 


f±(x,  t)  =  2rTi 
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where  q^(p)  are  tbe  roots  of  D(p,q(p))  *  0: 
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If  x  >  uQt,  the  first  term  in  Eq.  (10)  does  not  contribute,  since 
the  contour  may  be  closed  at  infinity  in  the  right  half-plane,  where  it  en*- 
closes  no  singularities,  (This  result  is  independent  of  whether  the  have 
only  simple  poles.)  The  expression  for  f^(x,t)  is  then  independent  of  the 
presence  of  the  boundary  and  consists,  as  in  the  unbounded  case,  of  plane  wave 
terms  proportional  to  exp  (p  (iV.)t  +  O.x)  where  p  (q) ,  l  -  are  the 

roots  of  D(p(q),q)  =  0. 

If  x  <  uQt,  the  first  term  of  Eq.  (10)  has  singularities  at:  (a) 

singularities  of  q  (p) ,  (b)  singularities  of 

X 

(a),  and  (c)  singularities  of  G^(p) ,  which  we  assume  to  be  simple  poles  at 
p  =  0j.,  i  -  1,.,,.  The  contour  can,  in  this  case,  be  closed  in  the  left  half¬ 
plane,  and  it  can  be  shown  that  the  singularities  (b)  exactly  cancel  the 
second  term  in  Eq.  (10).  This  cancellation  eliminates  the  plane  wave  terms 
which  appear  for  the  infinite  plasma.  In  their  ;  ace  the  singularities  (a), 

2  2  k 

which  include  the  essential  singularities  at  p  *  +  i  [uy  -  J* ,  give 

the  more  complicated  dependence  due  to  the  presence  of  the  boundary.  The 
singularities  (c)  give  rise  to  plane  wave  terras  proportional  to  export  + 

To  illustrate  the  long  time  behavior  of  an  initial  perturbation, 
consider  the  case:  np(x,0)  =  sinkx  and  £(x,0)  =  -N^ujp  k  (1-coskx) .  All 

other  boundary  and  initial  conditions  are  taken  to  be  zero.  The  integral 
(10)  around  the  essential  singularity  can  be  evaluated  by  the  saddle  point 
method  in  the  limit  of  large  t.  It  can  be  shown  that  this  method  is  valid 
in  the  present  case,  and  that  only  two  of  the  eight  saddle  points  which  occur 


H^(q^(p))  apart  from  those  in 


contribute  to  the  asymptotic  behavior.  The  result  is: 


S(x,t)  *  NQkuo2  (6tt)  ^  (uip/wB)2  T  5^6 

X  exp  v3  gX2/3  -  (%v/u>p)  t] 

X  sin  [0“3t  “  bX2/3r1/3  -  0  -  ^  tt]  (12) 

2  »i/5  * 

where  p  *  [1  -  (%v/u)p)  ]  ,  X  =  u)gX/uo>  t  “  wpt,  and  sin  #  *  This 

2/3  -2/3  2/3  1/3 

result  is  valid  provided  x  «  u  t  and  Px  T  ■«  1  «  px  T  It  should 

c 

be  noted  that  the  exponential  factor,  which  dominates  the  form  of  6(x,t),  is 

independent  of  the  initial  disturbance. 

The  above  results  yield  the  following  picture  for  the  development 

of  an  initial  disturbance.  If  the  disturbance  is  unstable  according  to  the 

usual  theory  for  spatially  unbounded  systems,  it  grows  as  e^t  only  until 

t  *  x/u  .  For  t  »  x/u  the  amplitude  goes  as  in  Eq.  (12).  It  reaches  a 
o  o 

maximum  at  tQ  *  P/y)3^2  (u^x/u^  and  finally  decays  as  e  ^ut.  If 

the  collision  frequency  is  large,  t  does  not  occur  in  the  asymptotic  region 

-i  f* 

t  >>  x/u  ,  so  the  behavior  is  always  e  in  this  region.  On  the  other 
o 

hand  if  v  =  0,  the  disturbance  grows  indefinitely,  and  one  has  essentially 

an  absolute  rather  than  a  convective  instability.  This  is  due  to  the  fact 

that  there  are  characteristics  of  Eqs.  (1-5)  which  are  parallel  to  the  t 
8 

axis.  The  presence  of  collisions  thereby  Insures  the  validity  of  the 
linear  approximation  for  all  times,  provided  x  is  not  too  large. 

In  experimental  measurements  what  is  observed  is  the  result  of 
many  disturbances  originating  at  different  prior  times.  Since  each  distur¬ 
bance  ha3  a  limited  lifetime  (for  any  fixed  x  and  v  #  0) ,  the  cumulative 
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effect  is  a  stationary  state  which  varies  with  x.  In  the  present  situation 
the  lifetime  depends  critically  on  the  collision  frequency,  so  one  expects 
the  level  of  fluctuations  in  the  stationary  state  to  also  be  strongly  in- 
fluc  :ed  by  v.  While  a  rough  estimate  of  this  stationary  state  can  be 
obtained  from  the  above  results,  a  more  accurate  analysis  is  desirable.  Such 
an  analysis  is  presently  being  performed. 
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